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Abstract 


The feasibility of using forebody tangential blowing to control the roll-yaw motion of a wind 
tunnel model is experimentally demonstrated. An unsteady model of the aerodynamics is 
developed based on the fundamental physics of the flow. Data from dynamic experiments is used 
to validate the aerodynamic model. A unique apparatus is designed and built that allows the wind 
tunnel model two degrees of freedom, roll and yaw. Dynamic experiments conducted at 45 
degrees angle of attack reveal the system to be unstable. The natural motion is divergent. The 
aerodynamic model is incorporated into the equations of motion of the system and used for the 
design of closed loop control laws that make the system stable. These laws are proven through 
dynamic experiments in the wind tunnel using blowing as the only actuator. It is shown that 
asymmetric blowing is a highly non-linear effector that can be linearized by superimposing 
symmetric blowing. The effects of forebody tangential blowing and roll and yaw angles on the 
flow structure are determined through flow visualization experiments. The transient response of 
roll and yaw moments to a step input blowing are determined. Differences on the roll and yaw 
moment dependence on blowing are explained based on the physics of the phenomena. 
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Introduction 


Advantages of flight at high angles of attack include increased maneuverability and increased lift 
during take-off and landing. At these flight conditions the aerodynamics of the vehicle are 
dominated by separated flow, vortex shedding and possibly vortex breakdown. These phenomena 
severely compromise the effectiveness of conventional control surfaces. Alternate means to 
control the vehicle at these flight regimes are therefore necessary. Several methods of active flow 
control have been proposed to provide this increased control'*. 

The present work investigates the augmentation of aircraft flight control system by the injection 
of a thin sheet of air tangentially to the forebody of the vehicle. This method, known as Forebody 
Tangential Blowing (FTB), is proposed as an effective means of increasing the controllability of 
aircraft at high angles of attack 5 7 . The idea is based on the fact that a small amount of air is 
sufficient to change the separation line on the forebody. As a consequence the strength and 
position of the vortices are altered causing a change on the aerodynamic loads. Celik 5 has shown 
that using this method side force, roll and yaw moments can be generated for angles of attack 
from 20 to 50 degrees. 

Experimental investigation is conducted to demonstrate the feasibility of using FTB to control 
the roll-yaw motion of a wind tunnel model. The model consists of a delta wing-body 
combination equipped with forebody slots through which blowing is applied. Experiments are 
conducted at a nominal incidence angle of 45 degrees 

A unique apparatus that allows the model two degrees of freedom, roll and yaw, is designed and 
built. The apparatus is used in dynamic experiments which show that the system is unstable, its 
natural motion divergent. External effects of the model support system are actively canceled. 

Flow visualization results revealed the vortical structure of the flow to be asymmetric even for 
the model at zero roll and yaw angles. The effects of blowing, roll and yaw angles on the flow 
structure are determined. Transient response of roll and yaw moments to step input blowing are 
characterized in terms of time constants. Differences on the roll and yaw moment behavior due to 
blowing are explained based on the physical mechanisms through which these loads are 
generated. 
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A model for the unsteady aerodynamic loads is formulated based on the basic physics of the 
flow. Parameters of the aerodynamic model are obtained from static and dynamic experiments. 
The aerodynamic model completes the equations of motion of the system which are used for the 
design of control laws using blowing as the only actuator. The closed loop control system is 
stable as is experimentally demonstrated in the wind tunnel. 


Experimental Apparatus 


Wind Tunnel and Model: 

The wind tunnel facilities of the Aeronautics and Astronautics Department at Stanford University 
are used for the experiments. It consists of a closed circuit low speed wind tunnel with 0.45m x 
0.45m test section. The maximum test section velocity is 50 m/s. 


As illustrated in Figure 1 , the model consists of a sharp leading edge delta wing with a 70 degree 
sweep angle and a cone-cylinder fuselage. Slots through which blowing is applied are present on 
both sides of the conical forebody. Air is provided to the forebody plena through flexible tubing 
that enters the model through the rear end of the fuselage. Tests are conducted at 20m/s and a 
nominal incidence angle of 45 degrees. The Reynolds number based on the wing root chord is 
260,000. Tunnel blockage at the nominal configuration is 7%. Figure 2 shows a side view of the 
test section with the model and the supporting mechanism. 



Figure 1 : Wind Tunnel Model and Detail of Blowing Slots 
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Model Support System^ 


A unique support system is designed and built to implement two degrees of freedom in the 
model 10 . The objective is to approximate the lateral-directional dynamics of an aircraft. Of 
particular interest is the roll-yaw coupling at high angle of attack. The apparatus can be divided 
into two subsystems: The first one implements the roll degree of freedom, <j>, and consists of a 
shaft mounted on bearings. The wind tunnel model is attached to the roll shaft allowing the 
model to rotate about its longitudinal axis. The roll subsystem is mounted on a mechanical arm 
that can rotate about an axis perpendicular to the models longitudinal axis (Figure 2). The 
mechanical arm implements the second degree of freedom, y, and is called the ‘yaw’ subsystem. 
For small roll angles, y equals the yaw rate. This approximation can be represented by relating 
<j> and y to the roll, pitch and yaw rates, p, q and r respectively. 


p = <j> q = ysin<)> r = ycos<|> 


( 1 ) 


Mechanical constraints limit the degrees of freedom as follows: l<|>l < 105° and lyl < 30°. 

An important aspect in the design of the experimental apparatus is that the dynamic properties of 
the support system should not dominate the dynamic response of the model. Experiments 
indicated that the friction in the bearings and the effect of the pressurized tubing are small when 
compared with the aerodynamic loads acting on the model*. 



Figure 2: Side View of Test Section and Two Degrees of Freedom Model Support System 
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Simulations and preliminary tests indicate that the inertia of the apparatus and the gravity 
restoring moment have a large effect on the overall system dynamics. Therefore an electro- 
mechanical system to cancel these effects has been designed and implemented 10 . The concept of 
active cancellation consists of providing a means to apply a torque that cancels the undesired 
external effects. For the current system a brushless motor is used to provide a torque which is 
computed based on measurements of the angle y and its angular acceleration. The same idea can 
be expressed in terms of the dynamic equations of the system: For the model wing at level, 4>=0, 
the equation of motion for the system can be written as: 

(I„ + UJ y - M n*.„ + (2) 

Where M represents the moment about the y-axis. 

For an ideal system: 


^model y ^«erod)lMBBC 

Therefore the torque that needs to be applied by the motor, M . is given by: 

M =1 v - M -M -M = I .7 - M (4) 

1VA motor 1 «nnl < AVA fnciian 1 A tnbu* gram? tnnl • gravity V 7 

The approximation is justified by the fact that for the current system, the inertia and gravity 
terms are dominant 10 . 

The inertia of the support, I^, is known, therefore by measuring the angular acceleration, the first 
term on the right hand side of the previous equation is obtained. The moment due to gravity is 
given by: 


M-* = — k^siny 


(5) 
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The constant kj, is a known quantity given by the product of the mass of the support by the 
distance of its center of gravity to the y-axis. The angle y is measured directly. In this way the 
torque that the motor should apply at any moment can be computed from these measurements. 

Figure 3 shows, in block diagram form, the implementation of the active cancellation loop. A 
torque sensor is designed and built to satisfy the specific requirements of this application. This 
sensor measures the torque that the motor applies to the system. A closed loop torque control 
logic is designed and implemented to provide fast torque command capability necessary for the 
cancellation of external effects. Two high precision linear accelerometers are used to provide 
measurement of the angular acceleration and a low friction precision potentiometer provides the 
measurement for the angle y. These two signals are fed into a computer that calculates the 
necessary torque command and sends it to the inner torque control loop that drives the motor. 



Figure 3: Block Diagram of Active Cancellation Loop. 


Air Injection System: 

Air can be injected independently through the slots located on both sides of the conical portion of 
the forebody. The amount of injected air is quantified by the jet momentum coefficient, which is 
defined as: 


C 




m jVj 


q-S 


ref 


( 6 ) 


Specially designed flowmeters are used to measure the mass flow rate from which Qi is 
calculated. Servo valves are used to vary the amount of injected air. The signal from the 
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flowmeters is read into a computer that implements a closed loop logic for C M control by 
commanding the servo valves. 


Instrumentation and Data Acquisition: 

The two degrees of freedom, (j) and y are measured by low-friction precision potentiometers. A 
six-component force-torque sensor connects the model to the roll shaft and is used to provide 
both static and dynamic measurements of the aerodynamic loads. Flowmeters measure the 
amount of blowing through each plenum. Two linear accelerometers are mounted to the y-axle. 
Their signals are combined to measure angular acceleration. A torque sensor connects the 
brushless motor to the y-axle and therefore provides a measurement of the torque that the motor 
applies to the system. 

For the flow visualization experiments, an argon-ion laser and an optical system are used to 
generate a laser sheet that is perpendicular to the model longitudinal axis. The optics is mounted 
on traversing system located on top of the test section allowing the laser sheet to be moved over 
the full length of the model. This capability is used in performing axial scans starting from the 
forebody and moving downstream to characterize the development of the flow structure. A 
smoke generator located upstream of the model is used to seed the flow. A video camera is 
located outside the test section aligned with the model longitudinal axis. The camera is used to 
record the results of the flow visualization and also the motion of the model during dynamic 
experiments. 

Three micro computers equipped with data acquisition boards are used in the experiments. One 
computer is dedicated to the active cancellation loop. A second computer is used to implement 
the closed loop control of the vehicle, i.e. to control the amount of air injected in each plenum. A 
third computer is used for data acquisition. 
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Experimental Results 


Flow Visualization: 

These experiments reveal the basic structure of the flow. Although four main vortices are 
expected, two from the forebody and two from the wing leading edges, experiments demonstrate 
that in general only three separate vortical structures can be clearly identified even for a 
symmetric condition in which <t>=y=0 and no blowing is applied (Figure 4b). 

By performing axial scans with the laser sheet it is observed that the asymmetry starts early on 
the forebody, i.e. close to the tip of the cone, and scales up over the entire forebody 10 . As a result 
of the asymmetry one vortex will be close to the fuselage and the other will be displaced and 



The effect of asymmetric blowing, i.e. blowing applied from one side only, is mainly to increase 
the asymmetry or invert it depending on which side the blowing is applied. Blowing moves the 
separation line on the forebody and can cause a change in the amount of vorticity that is shed. As 
a consequence the strength and positions of the vortices are affected by blowing. Experiments 
have shown that there is a finite amount of blowing that needs to be applied to invert the 
asymmetry of the flow, for the current model configuration at <J>=y=0 this value is C u = 0.0045. 

The application of symmetric blowing has the effect of changing the flow structure to a more 
symmetric one (Figure 5). For high values of symmetric blowing the flow can be considered 
attach ed on the forebody and its structure is very symmetric even on stations where the wing is 
present. 



(5a) Station 1 - Forebody - L/D=3 (5b) Station 2 - Wing-Body - L/D=5 

Figure 5: Smoke Flow Visualization Results - Symmetric Blowing, <|>=y=0, CV=.0075 

Experiments show that the asymmetric structure can change significantly by a change in the roll 
angle. The flow structure is not as sensitive to a change in y. 

Static Aerodynamic Loads: 

Static measurements for the roll and yaw moment as a function of <t>, y and are presented in 
Figures 6 through 8. A convention is adopted that right side, i.e. starboard blowing is positive 
and left side, i.e. port side blowing negative. In Figure 6 the effect of the roll angle on C, and C n 
is shown for y=0 and various C^. The C, curve for C,i=0 presents a change in slope for 4e-15°. 
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For Cji=.02 a change occurs at 4>=-5°. Also for C^=0 the C B curve presents a large change in slope 
for 4 k- 15°. The fact that these changes are not symmetric, i.e. they only occur for <(><0, indicates 
that they are caused by geometrical imperfections on the tip of the conical forebody. 




Figure 6: Roll Angle Effect on Roll and Yaw Moment Coefficients, C, and C„ , for y=0. 


Figure 7 presents curves for C, and C„ versus y for <|>=0 and various C M . Comparing Figures 6 and 
7 shows that the slopes of the roll and yaw moment curves are not as sensitive to y as they are to 
<|>. This is in agreement with the flow visualization experiments which indicate the flow structure 
to be less dependent on y than on <|>. 




Figure 7: Effect of y Angle on Roll and Yaw Moment Coefficients, C, and C n , for <J>=0. 

The effect of asymmetric blowing on C, and C B is shown in Figure 8 for y=0 and various <j). In 
this case blowing is applied either on the right or left side. As seen the roll moment varies 
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abruptly for IC M I<0.01 and the yaw moment presents a large variation for 0.01<Cn<0. Degani 9 has 
shown that small amount of blowing has an effect similar to the geometric imperfections near the 
tip of the forebody and can cause flow instabilities that lead to asymmetry in the flow. This is the 
most plausible explanation for the roll moment behavior shown in Figure 8. For IC^kO.Ol the roll 
moment is more sensitive to changes in than the yaw moment. This fact can be explained by 
the different physical mechanisms through which each of these loads are generated. For IC^kO.Ol 
the contributions of the direct jet momentum to C, and C„ are negligible. In this case, it can be 
considered that the roll moment is only generated at stations where the wings are present. 
Because the fuselage is of circular cross-section. Therefore the roll moment is determined by the 
flow over the wings which might be subject to large instabilities as vortex breakdown. The yaw 
moment is determined by the pressure distribution over the fuselage. The wings do not contribute 
to C„ because they are thin and only offer area perpendicular to the yaw axis. The flow over the 
forebody has a greater effect on the yaw moment due to its distance to the center of mass of the 



that occur in the asymmetric blowing case. The symmetry observed in these plots agrees with the 
symmetric flow structure obtained from the flow visualization experiments (Figure 5). 
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Figure 9: Effect of Symmetric Plus Asymmetric Blowing on Roll and Yaw Moment 
Coefficients for C^^O.Ol and y=0. 


Dynamic Experiments; 

The results for the static roll and yaw moment show that these moments are not zero for <t>=Y=0 
and 0^=0. Also the positive slope of the curve C o versus y indicates that the system is statically 
unstable at this condition. An experiment is performed to determine the dynamic characteristics 
of the system since those cannot be inferred from the static data alone: With no blowing applied 
the model is released from a certain initial condition. This represents the natural motion of the 
system and will ultimately determine if the system is stable or unstable. Figure 10 shows the time 
histories for <|> and y when the model is released from (Je'jsO, also shown are the results from 
simulati ons using the aerodynamic model described in the following section. As seen, the system 
is unstable. The motion is divergent and is stopped when the system approaches the mechanical 
limit of y. 
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Figure 10: Natural Motion of the Two Degrees of Freedom System. Experiment and Simulation. 


Equations of Motion 

In order to study the dynamic characteristics of the system it is necessary obtain its equations of 
motion. For the two degrees of freedom system those are: 

I M <|> + sin<|>cos<|>(I M -I M )y 2 + cos<(>l M y = M, 

x z y xz 

(1 A +I M sin 2 ()) + I Mi cos 2 (|>)y+ 2cos<(»sin<KI My -I Mz )4y+ (<^cos<t»-4> 2 sin<|))I M ^ = M 2 (7) 

For the current model configuration, where a vertical stabilizer is not present the product of 
inertia is zero. 


M, is the moment about the model longitudinal axis,. M 2 is the moment about the y-axis. M, and 
M 2 are given by: 


M, = M, A + M, t + M, P 

M 2 = M 2 a + M 2 t + M/ + M 2 ° + M 2 m (8) 


Where the superscripts indicates the origin of the moments: 
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A = aerodynamics 
T = air supply tubing 
F = friction of bearings and potentiometer 
G = gravity restoring moment 
M = motor 

For l<|>l < 40° the moment caused by to the air supply tubing on the first of equations (8) is 
negligible 10 . The torque applied by the motor is given by: 

M 2 m = - M 2 t - M 2 ° + I A y (9) 

The moments caused by friction of the bearings and potentiometers can be written as: 

M, P = - C F <j> M/ = -D P y (10) 

C p and D p are determined experimentally. 

Substituting expressions (9) and (10) in equation (8): 

M, =M 1 A -C F <(> 

M 2 = M 2 A -D F y + I A Y (11) 

Expressions for the aerodynamic moments M A and M 2 A are necessary to complete the above 
equations. Wong 3 developed an aerodynamic model for a delta wing undergoing roll oscillations 
that assumed that the dynamic loads could be approximated by lagged static loads and a pre- 
spedfied function of the roll rate. This basic idea is extended by including damping effects 
proportional to roll and yaw rates, cross-coupling terms in roll and yaw and apparent mass effects 
due to angular acceleration. The lag in the static loads is justified by comparing static and 
dynamic flow visualization results which clearly show that the vortex dynamic position lags with 
respect to the static one 7 . It is known that the strength of the shedding vortices is also affected by 
the motion of the vehicle. The current approach lumps position and strength effects by lagging 
the loads to represent their combined effect. Figure 11 shows the structure of the unsteady 
aerodynamic model in block diagram form. M, AS and Mj AS represent static roll and yaw moments. 
The unsteady effects are represented by the time constants x, and x 2 and functions f, and f 2 . 
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Figure 11: Structure of the Aerodynamic Model. M, AS and M^ s represent static 

roll and yaw moments. 

Functions f, and f 2 are assumed linear in the angular rates and accelerations: 


f i =C.<j> + C t y+C~$ 

f 2 = E>.^ + D t y + D ? y (12) 

Given the structure shown in Figure 1 1 and the expressions for f, and f 2 , the equations for the 
aerodynamic moments M, A and M 2 A are: 

Mj A = X, + C, +C^ + C t y + c^ 

m 2 A = *-2 + C 2 + d+ 9 + D fy + d?y (13) 


X, and X 2 are the lagged static roll and yaw moments for Qt=0, and are given by the following 
equations: 


TjX, +Xj =M^ 

't 2 X 2 +X 2 =Mf (14) 

C, and C 2 represent the effect of blowing on the roll and yaw moment and are given by: 
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^Ci+C,=Fi(QO 

^C 2 +C 2 =f 2 (Q0 


( 15 ) 


Where x$ and Zy are time constants that characterize the roll and yaw moment response to a 
variation in blowing, AC^. They include the effect of valve and plenum dynamics as well as the 
time it takes for the change in the flow structure to affect the moments. Figure 12 shows the 
response of C, and C„ for a step input command in Qi. Also shown are the results of the 
parameter fitting using equations (15). 




Figure 12: C, and C a response to Step Input Command in C^. AC, and AC n are used to 
indicate variation from the initial value. 


Substituting equations (1 1) through (15) into (7) the equations of motion for the two degrees of 
freedom system are obtained. Dynamic stability derivatives and time constants x, and x 2 are 
determined using a minimum least squares fit to the time histories of <)> and y for a set of dynamic 
experiments. Time constants z t and Zy are determined from the roll and yaw moment response to 
a step input command in Cji (Figure 12). 

The natural motion of the system is shown in Figure 10. Experimental results are compared with 
simulation using the aerodynamic model developed. The simulation agrees well with the 
measured response. 
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Closed Loop Control 


To design a control law that stabilizes the naturally unstable system the equations of motion are 
linearized as follows: For small static equilibrium roll and yaw angles, i.e. <!> E and y E « 1 , 4> and y 
are redefined as <t>-<|) E and y-y E and equations (7) and (1 1) are written as: 

I m J + d fT = m£ (16) 

About the static equilibrium position, M AS and M 2 A5 can be expressed as: 

Mj = C + 4> + C Y y 

= D^<|> + D T y (17) 

C* , Cy , D, and D T are the static stability derivatives obtained from the curves for the roll and 
yaw moments versus <J> and y. 

Equations (15) represent the effect of blowing. For the asymmetric blowing case and IQJ <0.01 
F, and F 2 are highly non-linear functions of , as seen in Figure 8. Furthermore operation in 
this region is prone to generating non-robust control laws because C, and C n are very sensitive to 
small variations in blowing. Given the above reasons the low blowing intensities are avoided by 
employing the following control strategy: A minimum amount of blowing other than zero, , 
is chosen and additional blowing AC 4 is added to that value, i.e.: 

Cn = Cy, + ACn (18) 

The curves for C, and C c versus AC H have the general form shown in Figure 13. 
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Figure 13: Characteristic of C, and C D versus AC^ 

A describing function approach is used to determine the equivalent gain of the curves C, and C„ 
versus AQ. The actual gain, N(A), depends on the amplitude of the input. A, and is given by: 


4D 

N(A) = + m 

kA 


(19) 


D and m are defined in Figure 13. An average amplitude is selected for AC^ and used to calculate 
C, and D b the equivalent gains for roll and yaw moment respectively. Expressions for F, and F 2 
are: 


F, = C B AQi. F 2 = D b ACji. (20) 

For the symmetric blowing case there is no need to apply the describing function approach 
because the C, and C„ variation with AC^.is fairly linear (Figure 9). 

The linearized equations of motion are: 

(ffci^ ~ C^)4» = (C^ — Cp)<j> + + Ci 

(I Mz “ D y)V = ( D ^ - D F)Y + D 4 < i > + ^2 + C2 

T,x, +X, =c^ + c 7 y 

^2^2 + ^2 = + 
t*Ci = — Ci +C B ACn 

T t C 2 = -C 2 + D B AC »* < 21 > 
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These can be written in the form: 


x(t) = Ax(t) + Bu(t) (22) 

Where x is the state vector and u the control variable: 

x(t) = [<t> 7 <)) 7 A,j A-2 Ci Cz] T 

u(t) = ACjx (23) 

A control logic is designed using the Linear Quadratic Regulator (LQR) method with weights on 
<(> , y, <(>, 7 and AC H . The result is a gain matrix, -K, that is multiplied by the state vector x(t) to 

generate the required control, i.e., 


u(t) = - Kx(t) (24) 

This control law requires knowledge of the state vector. <t> and y are measured directly. The other 
state variables are obtained from these measurements and the use of an estimator. 

The performance of the closed loop system is shown in Figure 14. The plots show data obtained 
during a real time closed loop control experiment. In this case the model was released from 
<|>=38 0 and ^14°. It is seen that the logic makes the system stable and regulates <j) and 7 to close 
to zero. The third plot shows the control effort, C^. Two curves are shown: Cn>0 for right side 
blowing and C^cO for left side blowing. 

Similar results are obtained using symmetric and asymmetric blowing. The disadvantage of the 
symmetric blowing case is a larger use of air. 
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Figure 14: Closed Loop Control of the Two Degrees of Freedom System. Response 
to Initial Condition. Asymmetric Blowing. 


Conclusions 

The use forebody tangential blowing (FTB) to stabilize the roll-yaw motion of a delta wing-body 
model is experimentally demonstrated in the wind tunnel. 

An aerodynamic model that is suitable for controls is developed based on: Static measurements 
of the aerodynamic loads and basic physical representation of the main dynamic effects. The 
model is validated through dynamic experiments and used in the design of closed loop control 
laws. The control logic stabilizes the system using blowing as the only actuator. It is shown that 
asymmetric blowing is a highly non-linear effector that can be linearized by superimposing 
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symmetric blowing. The transient response of roll and yaw moments to a step input blowing are 
determined. 

Dynamic experiments are conducted using a unique apparatus that allows a wind tunnel model 
two degrees of freedom, roll and yaw. These experiments show that at 45 degrees angle of attack 
the natural system is unstable presenting a divergent motion. 

The flow structure over the wing-body combination at 45 degrees angle of attack is asymmetric. 
As determined from flow visualization experiments. The coupling between forebody vortices and 
wing vortices is strong and an asymmetry that starts on the forebody will determine the structure 
of the flow downstream. At sections where the wing is present three main vortical structures are 
discernible. Asymmetric FTB increases the flow asymmetry or inverts it depending on which 
side of the model blowing is applied. The asymmetry can also be inverted by a change in roll 
angle. The flow structure is not as sensitive to changes in yaw angle. Differences on the roll and 
yaw moment dependence on blowing are explained based on the different mechanisms through 
which they are generated. 
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